Fisher Waves in the Diffusion-Limited Coalescence Process A + A ^ A 
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Fisher waves have been studied recently in the specific case of diffusion-limited reversible coa- 
lescence, A + A ^ A, on the line. An exact analysis of the particles concentration showed that 
waves propagate from a stable region to an unstable region at constant speed, just as in Fisher's 
"mean-field" theory; but also that the wave front fails to retain its initial shape and instead it 
broadens with time. Our present analysis encompasses the full hierarchy of multiple-point density 
correlation functions, and thus it provides a complete exact description of the same system. We 
find that as the wave propagates, the particles in the stable phase remain distributed exactly as in 
their initial (equilibrium) state. On the other hand, the leading particle — the one at the edge of the 
wave — advances as a biased random walk, rather than simply linearly with time. Thus the shape of 
the wave remains actually constant, but it is the "noisy" propagation of the wave's edge that causes 
its apparent broadening. 
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Fisher waves ||T|,|| are the best known paradigm of the 
invasion of an unstable phase by a stable phase. Fisher's 
theory describes the kinetics of particles concentration 
p{xi, X2, ■ ■ • , Xd, t) in a reaction process (taking place in 
d-dimensions) at the level of a reaction-diffusion equa- 
tion |,|: 



kip- k2p 



(1) 



Here D is the effective diffusion constant of the particles, 
and ki and ^2 are rates of particle generation and par- 
ticle death, respectively. Eq. (|^) admits two stationary 
solutions: 



and 



0, 



P = ki/k2 



(2) 



(3) 



but only the latter is stable. Suppose that the system is 
prepared at the initial state: p = ki/k2 for xi < 0, and 
p = for xi > 0. Then the stable phase {p = fci/fc2) 
invades the unstable phase (p = 0) in the form of a wave, 
p{xi, . . . , Xd, t) — f{xi — ct), which travels unchanged at 
constant speed c > Cmin — 2y/kiD. The minimal speed 
is realized for sufficiently sharp initial interfaces j^] (such 
as in our case). 

Fisher's equation might be regarded as a "mean-field" 
approximation to the kinetics of diffusion-limited coales- 
cence, A + A ^ A: The reaction terms of Eq. (|]) are 
in the form of the mass action rate equation appropriate 
for systems in local equilibrium, as might be expected for 
reaction-limited kinetics. Recently, however, diffusion- 
limited reversible coalescence was analyzed exactly in one 



dimension Q . It was found that, just as in the Fisher the- 
ory, waves propagate at a constant speed from a stable to 
an unstable phase, but that the width of the wave front 
broadens with time as w ~ ^/t. The broadening of the 
wave front is of special interest since it represents the ef- 
fects of the internal noise in the system — the noise which 
Fisher's equation fails to model. 

In ref Q| exact expressions were derived for the particle 
concentration, following the method of Intcr-Particle Dis- 
tribution Functions (IPDF). In this letter, we exploit the 
same method to derive the full distribution of particles, 
as represented by the infinite hierarchy of n-point density 
correlation functions. (The particle concentration corre- 
sponds to the special case of n = 1.) We find that as the 
wave propagates, the particles in the stable phase remain 
distributed exactly as in their initial (equilibrium) state. 
On the other hand, the leading particle — the one at the 
edge of the wave — advances as a biased random walk, 
rather than simply linearly with time. In this view, the 
shape of the wave remains actually constant in any par- 
ticular realization of the process. The apparent broaden- 
ing is the result of the fluctuations in the location of the 
wave's edge among different realizations (Fig. 1). 

The coalescence model l^-Q] is defined on a one- 
dimensional lattice of lattice spacing a. Each site is in 
one of two states: occupied by a particle A, or empty. 
Particles hop randomly into nearest neighbor sites, at 
rate D/a^. A particle may give birth to an additional 
particle, into a nearest neighbor site, at rate v/a (on ei- 
ther side of the particle) ||] . If hopping or birth occurs 
into a site which is already occupied, the target site re- 
mains occupied. The last rule means that coalescence, 
A + A — > A, takes place immediately upon encounter 
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of any two particles. Thus, together with hopping and 
birth, the system models the diffusion-limited reaction 
process A + A ^ A. The system's dynamical rules are 
illustrated in Fig. 2. 




(b) 

Figure 1. Fisher waves: (a) The ensemble average yields a 
wave traveling at speed c = v, with a front that broadens 
with time as w ~ V Dt. (b) Individual waves actually remain 
sharp throughout the motion, but their fronts travel in an er- 
ratic fashion, leading to the ensemble average in (a). Fisher's 
classical theory predicts a wave as in (a), but its front does 
not broaden with time. 
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(c) (d) 
Figure 2. Reaction rules: (a) diffusion; (b) birth; and co- 
alescence, (c) following diffusion, and (d) following a birth 
event. The dotted lines in (a) and (b) indicate alternative 
target sites. 

The IPDF method used for the exact analysis relies 
on the key concept of En,m{t) — the probability that sites 
n, n -|- 1, • • • , m are empty at time t. The probability that 
site n is occupied is 



Prob(site n is occupied) = 1 — E„ 



(4) 



The event that sites n through m are empty (prob. En^m) 
consists of two cases: site m -I- 1 is also empty (prob. 
En,m+i)i or it is occupied. Thus the probability that 
sites n through m are empty, but site m -|- 1 is occupied 
is En^m — En^rn+i- With this (and with a similar rule for 
when the particle is to the left of the empty segment) one 
can write down a rate equation_for the evolution of the 
empty interval probabilities 



9 (-^n,m — 1 Eji tn^ 



dEn^r. 

dt 
~-iE 



'n,m+l ) 



a- 
D 



2'(-E'ri+l,m 



En 



[(£'n,m — En^m+l) + {En,m ~ £'n-l,m)] • (5) 

0/ 

For example, the first term on the r.h.s. represents the 
event that sites n, . . . , m — I are empty and a particle at 
site m hops to m + 1, thus increasing -E„.m. The second 
term represents the decrease in En^m when a particle at 
m + l hops into the empty interval n, . . . ,m, etc. Eq. (|) 
is valid for m > n. The special case oi m = n yields the 
boundary condition 



En,n-1 — 1 



(6) 



The fact that the {En,m} represent probabilities implies 
the additional condition that En,m > 0. Finally, if the 
system is not empty then En^m as n —^ —oo and 
m —f oo. 

Rather than working with the discrete equations, it is 
simpler to pass to the continuum limit. We write x — na 
and y — ma, and replace -E'n.m(^) with E{x, y, t). Letting 
a 0, Eq. (||) becomes 

dE ^,92 92 



with the boundary conditions, 

E{x, x,t) = 1 
E{x,y,t)>0 
lim E{x, y,t) ^ 



(8) 
(9) 
(10) 



The concentration of particles is obtained using 
Eqs. (0) and (0), and passing to the continuum limit: 



dE{x,y,t) 
= \v=x 



(11) 



It can also be shown that the conditional joint proba- 
bility for having particles at x and y but none in between, 
is 



P2{x,y,t) = - 



d^E{x,y,t) 
dx dy 



(12) 



Given a particle at x, the probability that the next near- 
est particle to its right is at y, i.e., the "forward" IPDF, 
is 



Pix,y,t) = p{x,t) ^P2{x,y,t) 



(13) 
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Likewise, the "backward" IPD, the probabiUty that the 
next nearest particle to the left of a given particle at y is 
at X, is 

q{x, y, t) = t)-^P2{x, y, t) . (14) 
Eq. (^) admits two homogeneous stationary solutions: 
Eix,y) = l, (15) 

and 

E,^{x,y)=e~^^y--K (16) 

The first solution implies p — 0, while the second solution 
describes active equilibrium, with 



Pcq = V/D = 7 



(17) 



These states correspond, respectively, to the unstable and 
stable phases of Fisher's theory (Eqs. ^ and ||). It is 
important to notice that in the equilibrium phase the 
particles are distributed independently from each other. 
Indeed, suppose that the particles are so distributed, at 
a homogeneous concentration 7. Then, the probability 
that there are no particles in an infinitesimal interval of 
length A^, is (1 — 7A^). Since in equilibrium the particles 
are uncorrelated, the probability that a finite interval of 
length ^ = y — a; is empty, is (1 — 7A^)^/'^^. The equi- 
librium empty interval probability is then recovered by 
taking the limit Af — + 0. Notice also that in the equilib- 
rium state the IPDF (both forward or backward) is 



Pcq{x, y, t) = Qcqix, y, t) = 76 



(18) 



This can be derived from Eqs. (p^)-(p^), and (p^, as 
well as from the properties of uncorrelated, randomly dis- 
tributed particles. 

Suppose now that a system is prepared in the following 
initial state: at a; < the particles are distributed as in 
the equilibrium distribution, Eq. (^6|), and at x > there 
are no particles (the unstable steady state). The initial 
concentration profile is then 



p(x,0) = -[l-i/(.x)] 



v/D, X <Q 
0, 



a; > 



(19) 



where H{-) is the Heavyside step function. Doering et 
al. Q have shown that in this case the concentration 
profile evolves as 



p(x,t) = -^erfcf ^, ^ 



(20) 



where erfc(-) = 1 — erf(-) is the complementary error 
function That is, the initial step- wave propagates 

at speed v, while at the same time the wave front broad- 
ens as w ^ V Dt. 

Doering et al. have also considered the same initial 
condition as above, but when a particle is known for sure 
to be at the front's edge at x = 0: 



P(^,0) = -[1 



H{x)]+S{x) 



(21) 



In this case the wave front acquires an additional peak 
which propagates at the same speed and broadens in the 
same fashion as before: 



pix,t) 




y/AnDt 



{x - vtf 



(22) 



We now show that the distribution of particles in the 
coalescence system follows a simple and beautiful pat- 
tern: The particles in the stable phase — which are ini- 
tially distributed as in equilibrium (Eq. |l^) — remain dis- 
tributed in the same fashion throughout the lifetime of 
the wave. The noisiness of the system is manifested in 
the motion of the leading particle (the one at the wave's 
edge): rather than moving at constant speed c = v, \t 
performs a biased random walk with average drift veloc- 
ity V. Thus, the broadening of the wave front is the re- 
sult of averaging over the shifting positions of the leading 
particle in different realizations of the process (Fig. 1). 
For simplicity, we first show that such an interpretation 
is consistent with the exact solution of ref. [Q, and we 
defer the complete proof until later. 

Let the position z of the leading particle be given 
by the distribution probability density p{z,t). Let the 
other particles — those to the left of the leading particle — 
remain distributed as in equilibrium at all times. The 
probability that the interval {x, y) is empty at time t de- 
pends upon the location of the interval endpoints, x and 
y, with respect to that of the leading particle, z: 



E{x,y,t) 




z < X < y 
X < z < y 
X < y < z 



(23) 



Hence, taking into account the distribution of z, 
E{x,y,t)= p{z,t)dz + e-^^y"'n p{z,t)dz 



= P(x,t)-)-e-^(^--)[l-P(y,t)] , (24) 
where in the last equation we introduced the definition 

Piz,t) = [ p{z',t)dz' . (25) 



In terms of p, the concentration of particles (Eq. |ll[ ) be- 
comes 



p{x,t)^p{x,t) + ^[\-P{x,t)] 
Notice also that the forward IPDF is 



(26) 



p{x y t) = ^e-7fo-.)P(yito£^J^Ml (27) 
^ ^ p{x,t)+^[l-P{x,t)]^ ^''^ 
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while the backward IPDF remains exactly as in equilib- 
rium (Eq. |l8|) , consistent with the postulated distribution 
of particles. 

It is straightforward to verify that the proposed form 
of E(x,y,t), Eq. (p^), satisfies the boundary conditions 
(i)-®- Putting E{x,y,t) in Eq.(0) we sec that it is 
also satisfied, provided that 



d 



d 



P{z,t)^D-^Piz,t)~v—P{z,t). (28) 
at oz'^ oz 

From the definition of P, we have the boundary condi- 
tions: 



lim P{z, t)^0 , 
lim P{z,t) = 1 , 
P{z,t) > . 



(29) 
(30) 
(31) 



Finally, if a particle is known for sure to be present ini- 
tially at X = (Eq. |2^) , that particle is clearly the leading 
particle, and p{x, 0) = 5{x). Thus, 



P{z,0) = H{z) 



(32) 



The solution to Eq. (|28| ) which satisfies the boundary 
conditions (p9|)-(^ll) and the initial condition (^2|), is 



Piz,t) 



(33) 



Thus, using Eq . (|26| ), we recover the prediction of Doer- 
ing et al., Eq. (|22|). The probability density function for 
the position of the leading particle is particularly simple: 

p{x, t) ^ {AiiDt)-'^/^ exp[-(x - vtf/ADt] , (34) 

and identical to that of a Brownian particle character- 
ized by the diffusion coefhcient D and subject to a drift 
V. A moment's reflection shows why this is the case: The 
particle at the edge of the wave can step to the right or 
left with equal probabilities, at rate D/o?. This explains 
diffusion. The drift is a result of the birth mechanism. 
If the leading particle gives birth onto the site to its left, 
the edge does not move. However, if the particle gives 
birth onto the site to its right then the edge moves to the 
right (at rate v/a). Because of this left/right asymme- 
try the walk performed by the edge is biased. — What 
is surprising, is that the particles trailing the edge re- 
main distributed, on average, exactly as in their initial 
equilibrium distribution! 

Consider now the initial condition (|l9|), where the lead- 
ing particle is not necessarily at z = 0. Because the parti- 
cles are distributed as in equilibrium, the leading particle 
is at z < with probability p(z, 0) = 'ye^^ (and at z > 
with probability 0). This case may be regarded as a su- 
perposition of systems of the previous type, where the 
leading particle is initially surely at z. Indeed, the linear 
combination 



p{x,t) 



^e^-{_Lerfc(^^^^) 



1 



V47rDi 



exp[-i:^^^]}dz, (35) 



reproduces the known result of Eq. ( pCi| ) . 

Up to this point we have merely found a distribu- 
tion of particles which happens to explain the known re- 
sults for E and for the concentration of particles p. It 
is conceivable, however, that other distributions might 
accomplish the same feat. A complete description of 
a particle system requires knowledge of the full hier- 
archy of n-point density-density correlation functions 
Pn{xi,X2, . ■ . ,Xn,t): the joint probability to find parti- 
cles at xi,X2, ■ ■ ■ jXn at time t. The particle concentra- 
tion which was studied above corresponds to the special 
(limited) case of n = I . The n-point correlation functions 
can too be analyzed exactly through the IPDF method 
1^. We shall now describe the procedure and employ it 
for the conclusion of our proof. 

Let En{xi,yi,X2,y2, ■ ■ ■ ,Xn,yn,t) be the joint proba- 
bility density that the intervals {i = l,2,...,n) 
are empty at time t. The intervals are non-overlapping, 
and ordered: xi < yi < ■ ■ ■ < Xn < y-n- Then, the 
n-point correlation function is given by 



Pn ('^l 1 ■ ■ ■ ; ; ^) 



dyi--- dy, 



-En{xi,yi, . . . ,Xn,yn,t)\yi=Xi ■ (36) 



Doering |^ has shown that in our coalescence system the 
En satisfy the partial differential equation: 



d 

— En{xi,yi, . . .,Xn,yn,t) 
Ot 



D(— — 

dx\ dyl dxl dyl 



Er, 



d d d d ^ 

dxi dyi dxn dyn 

with the boundary conditions 

lim En{xi,yi,. . . ,Xn,yn,t) = 

Xi^^yi or Uiixi 

En-i{xi,yi, ....i^i.ii,.. .,x„,yn,t) , (38) 

and 

lim En{xi,yi, . . . ,Xn,yn;t) ^ 

yilxi+i or Xi+iij/i 

En-i{xi,yi, . . . , ^i, ■ ■ .,Xn,yn;t) , (39) 

which should be obeyed for all n > 1. Here, we use the 
notation that crossed out arguments (e.g. f^i) have been 
removed. Notice how the E„ are tied together in an hi- 
erarchical fashion through the boundary conditions ( [3^ ) 
and (|39|): one must know E„-i in order to compute En- 
At the root of the hierarchy, Ei = E is the simple empty 
interval probability that was studied above. 
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If our interpretation of the particles distribution is cor- 
rect, then, following a reasoning similar to that which led 
to Eq. (|2|) , we should have 

+ ■■■ + e-^««i--i)+-+(^---')>[P(x,+i,t) - P{y,,t)] 
+ --- + e-''«^i-^i)+-+(^"-^")>[l-P(2/„,i)] . (40) 

It is easy to confirm that these functions fulfill the bound- 
ary conditions (jS^) and (|39|). Eq. ( |37| ) is also satisfied, 
provided that P satisfies the same equation as before, 
Eq. (H). Using Eqs. (||), (||), and @, we find the 
n-point correlation function: 



.0 =7""V(a;„,0 



(41) 



Thus, the joint probability density for finding particles 
at xi < X2 < ■ ■ ■ < Xn depends only on the rightmost 
particle at Xn — it alone is represented by the overall con- 
centration profile. The particles ai Xi, . . . , Xn-i are dis- 
tributed randomly and independently of each other, with 
the equilibrium density 7. This is exactly what one would 
expect from the particle distribution that we have cham- 
pioned throughout this work, and it completes the proof 
to our claim. 

In summary, we have shown that in one dimension the 
broadening of the wave front of Fisher waves in the coa- 
lescence system is truly due to the spread of locations of 
the leading particle in different realizations of the process. 
This peculiar situation, where the particles remain dis- 
tributed as in the initial equilibrium configuration — and 
where all the changes are only reflected in the position of 
the leading particle — was recently found also in diffusion- 
limited coalescence in one dimension in the presence of a 
trap In fact, the proof in the present article may be 
carried over to that situation with only minor changes. 

Fisher waves in the coalescence system have been stud- 
ied numerically also in higher dimensions 1 11 1 . The goal 



an exponent which decreases with increasing d and be- 
comes zero (as in Fisher's mean-field description) above 
the critical dimension of rf = 3. The width of the wave 
front in ref. []lT| was measured from the concentration 
profile, as averaged over an ensemble of different realiza- 
tions. It would be interesting to repeat these simulations, 
but this time to subtract the ffuctuations in the position 
of the wave front between different runs, as suggested 
from the present findings. This could result in a decrease 
of the width exponent, and possibly also in a decrease of 
the critical dimension for the validity of the mean-field 
Fisher picture. 

I thank Paul Krapivsky for stimulating discussions. 



there was to explore the effect of fluctuations intrinsic 
to particle systems as a function of dimensionality. The 
wave front was found to broaden as a power of time, with 
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